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Abstract 

We present preliminary results for the measurement of the branching fraction of the decay 
XcqK^ from a sample of 75 million BB pairs collected by the BABAR detector at the PEP-II 
asymmetric-energy B Factory at SLAC. The XcO meson is reconstructed through its two-body 
decays to vr+vr- and K+K~ . We measure B{B^ XcoK^)xB{xcO ^ vr+vr-) = (1.08±0.35(stat) ± 
0.10(syst))xl0-6 and B{B^ XcoK^)xB{xcO ^ K+R-) = (1.48±0.44(stat) ±0.17(syst)) x IQ-^. 
Using the known values for the XcO decays branching fractions, we combined these results to obtain 
B{B^ XcoK^) = (2.4 ± 0.7) x lO"! 
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1 Introduction 



The XcO state corresponds, in spectroscopic notation, to the 1 Pq state of the charmonium spectrum 
(jPc: = 0++), whose mass and width are measured to be (3415.1±0.8) MeV/c^ and (14.9;2:3) ^eV, 
respectively [||. 

In the framework of the factorization approximation, the production rate for B XcqK decays 
vanishes when the assumption is made that this process occurs through the decay of a 6 quark 
into a color-singlet cc pair at short distances Q. It has been suggested |^, |^ that the "color- 
octet mechanism" could enhance such decays via the emission of a soft gluon. The branching 
fraction for B XcO ^ would then be expected to be comparable in magnitude to that of the 
"color-singlet-allowed" decay mode B — > Xci K. 

Belle recently reported the first evidence for the B^ x^qK^ decay On a sample of 
31.3 X l{f BB events, Belle measured B{B^ XcqK^) = [Q-^tW ± 1.1) x lO""^; the XcO meson 
was reconstructed through its decay to the vr+vr" mode. Previously CLEO had reported an upper 
limit of B{B^ XcoK^) < 4.8 x 10"^ at 90% C.L. |]. 

The results of a similar analysis performed on BABAR data are reported in this paper; both 
XcO vr+vr" and XcO — ^ K^K~ decay modes have been used for the measurement of B{B^ 
XcqK^); the branching fractions of XcO ^ tt+tt" and XcO K^K~ are (5.0 ± 0.7) X 10-3 and 
(5.9 lb 0.9) X 10^3 respectively ||l|. The choice of these two modes was motivated by their cleaner 
signature compared to other higher multiplicity hadronic modes showing larger branching fractions. 

2 The BaBAR detector and dataset 

The data used in this analysis were collected with the BABAR detector at the PEP-II e^e~ storage 
ring. The data sample we used corresponds to about 70 fb~^ of integrated luminosity collected 
at the T{4:S) resonance ( "on-resonance" ) between October 1999 and May 2002. The collider is 
operated in an asymmetric mode, that is with beams of unequal energies, so that the center of 
mass is "boosted" (with (3^ = 0.55) along the direction of the electron beam. 

The BABMi detector is fully described elsewhere 0. It consists of a charged particles tracking 
system, a Cherenkov detector (DIRC) for particle identification, an electromagnetic calorimeter and 
a system for muon and identification. The tracking system consists of a 5-layer, double-sided 
silicon vertex tracker and a 40-layer drift chamber (filled with a mixture of helium and isobutane) , 
both in a 1.5 T magnetic field supplied by a superconducting solenoidal magnet. The DIRC is 
a unique imaging Cherenkov detector relying on total internal reflection in the radiator. The 
electromagnetic calorimeter consists of 6580 CsI(Tl) crystals. The iron flux return is segmented 
and instrumented with resistive plate chambers for muon and identification. 

3 Analysis method 

A preliminary selection for hadronic events required the presence of at least three charged tracks, 
the ratio between the 2nd and 0th order Fox-Wolfram moments to be lower than 0.5 and the 
total energy of all the charged and neutral particles in the event to be greater than 4.5 GeV. 
Furthermore, at least one track identified as kaon was required to have a momentum greater than 
900MeV/c in the center of mass frame. 

B^ meson candidates were reconstructed by combining a track identified as a charged kaon 
with a XcO candidate. The kaon four-momentum was required to be loosely consistent with the 



8 



kinematics of the decay; the XcO nieson was reconstructed from a pair of charged tracks identified 
as both pions or both kaons. The typical efficiency for the kaon selectors used is in the range 70% 
to 90%, depending on momentum, while the typical probability for a pion to be misidentified as a 
kaon is below 5%. All these tracks were explicitly required to have polar angles lying in the region 
0.35 rad < 9 < 2.54 rad, to have at least 12 hits in the drift chamber and a transverse momentum 
with respect to the beam direction greater than 100 MeV/ c. In addition, tracks identified as decay 
products of Kg vr+vr^, r] tt^tt^tt^, A ^ pn ot tracks from 7 conversions were rejected. The 
XcO candidates were then defined as those for which the Tr~^TT^ {K^ K ) invariant mass lies in the 
region 3.32 GeV/c^ < m^^o < 3.50 GeV/c^. 

In order to reject the large background coming from continuum events, a Fisher discriminant 
has been used, built from a linear combination of eleven event shape or B kinematics related quan- 
tities. The coefficients were determined, separately for the two modes, by maximizing the achieved 
separation between signal and continuum background (both from samples of events simulated with 
Monte Carlo techniques). 

The selection of B candidates relied on the kinematic constraints given by the T(4S') initial 
state. Two largely uncorrelated variables were used: the beam-energy substituted mass, mss = 
\/^beam " |P*P) and AEev = E* - \J m\ + |p*|2, where El^^^ = ^/2 is the expected B energy 
in the T{AS) rest frame, E* and p* are respectively the energy and vector momentum of the B 
candidate in the T(4S') rest frame, and is the mass of the B^ meson. In cases where multiple B 
candidates were present in the same event, only the one with the smallest absolute value of Ai?ev 
was retained. 

The values of the cuts for the Fisher discriminant, rriEg and AE'ev were determined by an 
optimization procedure aimed at maximizing the value of S/^J S + B. The number S of signal 
candidates and B of background events surviving the cuts were estimated on samples of simulated 
events and on data from the AE'ev "sidebands" of the rriES vs. AE'ev plane, respectively. The 
sidebands were defined by 5.2GeV/c2 < < 5.3GeV/c2, 0.1 GeV < jAE^evI < 0.2 GeV. The 

signal and background samples were normalized to each other assuming the value measured by 
Belle for B{B^ XcoK^) and the world average for B{xcO tt+tt^) and B{xcO K^K^) Q. 

The signal region in the rriES vs. AE'ev plane was defined by —45 MeV < Ai^cv < 35 MeV, 
mES > 5.275 GeV/ for the XcO ^ Tr+vr" mode and by -70 MeV < l^E^^ < 60 MeV, mES > 
5.2735 GeV/c^ for the XcO K~^K~ mode. To account for a shift in AE'ev between Monte Carlo 
and data, observed on samples of 5 ^ D^tt events with a vr+vr^ final state, a correction of +7 
MeV has been applied to Ai?ev on data for both modes. Table [l| summarizes these cuts and the 
number of selected events in the signal box. 

Table 1: Signal box definition (including the region in iti^^q) and number of selected events for the 
two modes. 



Mode 


m^^o(GeV/c^) A^ev(MeV) mEs(GeV/c^) 




XcO vr+vr- 
Xco ^ K+K- 


3.32 ^ 3.50 -45 35 > 5.2750 
3.32 ^ 3.50 -70 60 > 5.2735 


83 
84 



In order to remove the contamination from B decays to D and D* modes, an explicit veto 
on several of these modes has been applied: B candidates were rejected if at least one of their 
decay products also contributed to the reconstruction of a i? in one of the veto modes, with 
|A£;ev| < 30 MeV and tues > 5.27 GeV/ c^. 
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Some residual contamination was still observed from modes with a meson decaying to K^n^ 
OT a 4> meson decaying to K'^K~, for — *■ XcO 7r~^7r~)K^ and — XcO K^K~)K^ , 
respectively. This background was suppressed by rejecting all events for which the invariant mass 
of the pair formed by the fast kaon with the oppositely charged pion (kaon) was compatible with 
the {(p) mass. 

The main source of non-combinatorial background remaining after the selection described comes 
from non resonant B decays with the same final state as the signal: 5=*= — K"^ tt+tt" and 5=*= — 
K^K~ . The fraction of such events surviving the set of cuts detailed above is 1.4 x 10~^ and 
3.5 X 10~^, respectively; these estimates however could not be used for a reliable evaluation of the 
expected contamination, since the associated branching fractions are not well known. These modes 
are expected to peak in mES and AE'ev (what will be referred to as "peaking background" in the 
following), while the distribution of m^^^^ is expected to be flat: this was used to separate their 
contribution from the signal by means of a fit to the data, as will be described in the following 
section. 

The background from mis-reconstructed XcO decays to other modes has been estimated to be 
negligible with respect to the other background sources for both the vr+vr" and the K~^K~ modes. 

The overall selection efficiency, estimated from Monte Carlo, is (31.2 zt 0.6 zt 1.6)% for the 
XcO 7r+7r~ mode and (24.1 it 0.5 ± 1.1)% for the XcO K'^K^ mode. In both cases, the first 
error is statistical and the second reflects the systematic uncertainties. These are mainly due to 
differences between data and Monte Carlo in the track reconstruction and particle identification 
efficiencies; a smaller contribution comes from differences in the efficiency of the selection cuts, in 
particular the cut on the Fisher discriminant. 

4 Extraction of the signal yield 

The extraction of the signal yield was complicated by the large fraction of background events (both 
combinatorial and peaking background) in the selected sample. B decays were separated from 
the combinatorial background by means of the two kinematic variables tties and AE'ev; the m^^^g 
distribution was used to separate the signal and the peaking background contributions. 

The number of signal events was extracted by fitting the m^^-^,, distribution for events in the signal 
box with a generic background component (combinatorial + peaking) and a signal component. An 
unbinned maximum likelihood fit was used. 

Background events were assumed to have a uniform m^^^ distribution in the narrow region 
selected for this variable. This was confirmed by looking at the m^^^ distribution for data events 
away from the (msS) Ai?ev) signal box, and for B^ —y 7r"'"7r" - B^ ^ K+K- Monte 
Carlo events. The signal component was modelled by a non relativistic Breit-Wigner function 
convoluted with a Gaussian function. In the fit the mass and the total width of the XcO were fixed 
to the expected values; the width of the Gaussian function was determined from Monte Carlo. The 
reliability of the fit was checked on a "cocktail" sample containing known amounts of combinatorial 
background, peaking background and signal events. 

The resulting signal yields are A^sig(XcO ^ vr+vr") = 25.4±8.2 and iVsig(XcO ^ K+K') = 27.0± 
8.0. The breakdown of the background contributions, as estimated by combining the above results 
with fits to the rriES distributions, is -/Vbkg(XcO — ^ 7r"'"7r~) = 26.1(combinatorial)+31.5(peaking) and 
A?bkg(XcO — *■ K'^K~) = 20.8(combinatorial)-|-36.2(peaking). These values refer to the whole region 
3.32 GeV/c2 < m^^^ < 3.5 GeV/c^. 
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Figure [T| shows the rriEs distribution for events in the AE'ev signal region and the fn^^^g dis- 
tribution for events in the signal box, for the two modes considered. Superimposed are the fit 
results. Comparison of the thes and m^co plots reveals the presence of a substantial component 
of non-XcO events with a peaking structure in mES in both cases. Furthermore, the enhancements 
seen in the m^^,, distribution appear shifted to lower values with respect to the expected XcO mass. 





Figure 1: Extraction of the signal yield for XcO ~^ it~^tt~ (top) and XcO K~^K~ (bottom). Left: 
distribution of the mES variable for events selected in the AE'ev signal region. Superimposed is the 
result of the fit to a Gaussian B signal component (solid line) + an Argus function for the continuum 
background (dotted line). Right: distribution of the recoil invariant mass m^^^^ for events in the 
signal box. Superimposed is the result of the fit to a Breit-Wigner function convoluted with a 
Gaussian function (solid line) + a flat background component (dotted line). 
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The distortion of the m-^^p distributions could originate from the presence of interference effects 
with non-resonant 3-body decays; this effect has not been taken into account in this analysis. As 
a check, we repeated the fit by allowing the XcO mass to float, and obtaining for it values shifted 
by about -6 MeV/c^ (for XcO vr'^'Tr") and -10 MeV/c^ (for XcO K^K~) with respect to the 
published PDG values. Results for the yield were compatible well within 1 a with those extracted 
with the mass fixed. 

As a way to estimate systematic effects related to the fitting procedure, an alternate method was 
employed. There, the mss and m-^^^ distributions were simultaneously fit for all events in the AE'ev 
signal region. In this case the likelihood function contained three terms: one for combinatorial 
background, one for the peaking background and one for the signal. The ttt-es distribution for 
the peaking background component was modelled with a Gaussian with the same resolution as 
the signal, while the combinatorial component was fit to the usual empirical Argus shape. Both 
backgrounds were assigned uniform distributions in m^^^Q. This fit yields A^sig(XcO vr"'"7r~) = 
23.4 lb 8.5 and Nsig{xcO K~^K^) = 29.7 it 8.5. The absolute difference in the signal yield with 
respect to the previous method is taken as a global systematic error. 

5 Branching fraction measurement 

The branching fractions are obtained as B = where e denotes the overall signal efficiency 

and Ng± is the total number of mesons produced in the data sample considered. 

Nq± is obtained from the measured number of BB pairs, N^-^ = (75.7 it 1.1) x 10^, assuming 
e(r(45) ^ B+B-) = 50%. 

The results for the product of the branching fractions for the two modes are reported in Table |2|. 
The systematic errors reported there combine the uncertainty on the efficiency evaluation with the 
error on and that associated with the fitting procedure. 

In order to extract the value of the branching fraction for B^ XcoK^, we used the values 
reported by the PDG [||] for B{xcO 7T~^7r~) and B{xcO — > K~^K~), to obtain: 

B{B^ XcoK^) = (2.15 ± 0.69(stat) ± 0.21(syst) ± 0.30(syst)) x 10"^ 

for XcO ^ vr"*"7r~ and 

B{B^ XcoK^) = (2.51 ± 0.74(stat) ± 0.28(syst) ± 0.38(syst)) x lO""^ 

for XcO ~^ K^K~; here the second systematic error accounts for the uncertainty on the XcO 
decay branching fractions. The two measurements were combined to obtain: 

B{B^ XcoK^) = (2.4 ± 0.7) x 10"^ 

When combining the errors, we treated the systematic uncertainties arising from efficiency 
corrections, from the fitting procedure and from the number of BB pairs determination as fully 
correlated, while statistical errors have been treated as uncorrelated, since no events belong simul- 
taneously to both the samples. To combine the errors coming from the measurement of the XcO 
branching fractions, we only used the values published by the BES Collaboration and treated 
the errors as fully correlated. 

The measurement obtained for B{B^ XcoK^) is of the same order of magnitude as that 
for B^ — > XciK^ ((6.5 it 1.1) x 10~^ from [l|), thus supporting the treatment of charmonium 
production in B decays with a QCD factorization approach including non-perturbative color-octet 
contributions. 
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Table 2: Branching fractions from the two XcO decay modes and their combined value. The meaning 
of the quoted errors is explained in the text. 



Mode 


B{B^ ^ XcoK^) X Bixco ^ /) (10-«) 


B{B^ ^ XcoK^) (10-4) 


XcO TT+vr" 
XcO ^ K+K- 


1.08 ± 0.35(stat) ± O.lO(syst) 
1.48 ± 0.44(stat) ± 0.17(syst) 


2.15 ± 0.69(stat) ± 0.21(syst) ± 0.30(syst) 
2.51 ± 0.74(stat) ± 0.28(syst) ± 0.38(syst) 


Combined 




2.4 ±0.7 
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